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ABSTRACT 

A logarithmic transform of the convergence field improves 'the information content', i.e., the overall pre- 
cision associated with the measurement of the amplitude of the convergence power spectrum by improving 
the covariance matrix properties. The translation of this improvement in the information content to that in 
cosmological parameters, such as those associated with dark energy, requires knowing the sensitivity of the 
log-transformed field to those cosmological parameters. In this paper we use A^-body simulations with ray 
tracing to generate convergence fields at multiple source redshifts as a function of cosmology. The gain in in- 
formation associated with the log-transformed field does lead to tighter constraints on dark energy parameters, 
but only if shape noise is neglected. The presence of shape noise quickly diminishes the advantage of the log 
mapping, more quickly than we would expect based on the information content. With or without shape noise, 
using a larger pixel size allows for a more efficient log-transformation. 

Subject headings: cosmology: theory - gravitational lensing large-scale structure methods: numerical 



1. INTRODUCTION 

The nature of the dark energy is one of the most intrigu- 
ing mysteries of the Universe. As a result, various large sky 
area surveys are being conducted and designed to statistically 
determine the properties of this energy component with high 
precision. We of course wish to extract the most informa- 
tion available from the data. Recently, it has been suggested 
that the two-point statistics of the logarithmically transformed 
nonlinear density or weak-lensing convergence field may 
contain more information than the c onventional nonlinear 
fields without t he transformation (e.g., iNevrinck et al.ll2009t 
ISeo et al]|20I Ih . In linear theory (i.e., at high redshift), the 
convergence (and density) field is Gaussian, which means that 
the two-point function contains all the information. Due to 
the structure growth, the field however becomes more nonlin- 
ear and non-Gaussian at low redshift. The cosmological in- 
formation in the two-point function therefore decreases with 
increasing nonlinearit y, the lost informati on moving to the 
higher order statistics (iTakada & Jainll2004l) . 

A logarithmic transform of the nonlinear mass/galaxy den- 
sity field or the weak-lensing convergence field, which makes 
the one-point distribution of the field r nore Gaussian (e.g ., 
iColes & JonesI [19911 iKavo et al.l l200Tt iTaruva et al.lllool . 
appears to produce a final field that alleviat es this problem 
by mimicking properties of a Gaussian field (INevrinck et al.l 
[2009; Seo et al. 2011; Yu et al. 2011). The two-point func- 
tion of the t ransformed field has a more diagonal covariance 
matrix (e.g.. INevrinck et ani2011al) and (therefore) increased 
information content (i.e., the overall precision associated with 
the measurement of the amplitude of the convergence/density 
power spectrum) to a level comparable to the Gaussian field. 
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ISeo et all (1201 ll) showed that a Taylor expansion of the loga- 
rithmic transformation that includes up to the bispectrum con- 
tribution captures most of the improvement on large scales, 
suggesting that the log-transform draws its extra information 
from higher order statistics. 

In these works, the benefits of the log-transform have been 
studied only when a single parameter - the amplitude of the 
power spectrum - is varied. This does not necessarily trans- 
late into an improvement in the measurements of other cos- 
mological parameters, such as dark energy parameters. If 
the power spectrum of the transformed field is less sensi- 
tive to, e.g., dark energy (i.e., smaller derivatives with re- 
spect to the parameters), or suffers more degeneracies be- 
tween parameters, the log-transform may not be as effec- 
tive as expected based on the improved information content. 
In order to test this, we need to understand the dependence 
of the log-transformed field as a function of cosmology. In 
ISeo et al . (2011), which we refer to as Paper I hereafter, we 
used a modified log-transform for the weak-lensing conver- 
gence field and showed the increased information content af- 
ter the log transform. Here we extend the previous work and 
investigate whether or not this improvement propagates to the 
determination of the dark en ergy parameters . A fi rst step 
in this direction was taken in iJoachimi et al.l (1201 ll) for the 
weak lensing field, wherein mapping the required derivatives 
with respect to cosmological parameters were computed ana- 
lytically. They showed that the Box-Cox transformation that 
encloses the logarithmic transformation, when opti mized, in- 
deed gi ves a tighter constraint on 0„i — erg. Also. iNeyrinckl 
(1201 Id) has recently shown tha t, using Coyote Universe sim- 
ulations dHeitmann et al.llTOlOl) , log-transformed density field 
gives tighter, but unmarginalized, constraints on cosmological 
parameters. 

Here, we adopt a more extensive numerical approach and 
derive the dependence of the transformed field on the cosmo- 
logical parameters directly from a very large set of A^-body 
simulations produced for various cosmologies. The total sim- 
ulated sky is 37000 square degrees. We produce mock weak- 
lensing convergence fields at three source redshifts, carry out 
a Fisher matrix analysis of tomography using the fully non- 
linear covariance matrix and derivatives, and derive marginal- 
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ized errors on 6-8 cosmological parameters, including dark 
energy parameters {i.e., fix, wq, and Wa)- As far as we know, 
this paper is the first to conduct the full Fisher matrix anal- 
ysis of weak-lensing power spectrum tomography using the 
numerical simulations to compute the derivatives and covari- 
ance matrices (but also see iYang et al.i201 U for peak statistics 
tomography). 

The paper is organized as following. In § |2] we explain the 
details of our A^-body simulations and Fisher matrix analysis. 
In § [3] we revisit the general properties of the log transform, 
such as the one-point probability distribution function (PDF), 
the information content, and the covariance matrix from the 
A^-body results using three source redshift bins. In § |4] we 
present the numerical derivatives before and after the logarith- 
mic mapping and the results of the full Fisher matrix analysis, 
deriving the dark energy figure-of-merit. We compare this re- 
sult with the Fisher matrix obtained from semi-analytic fits to 
the power spectrum. In § |5] we include shape noise; in § |6] 
we discuss the effect of the size of pixels; in §|7] we discuss 
the analytic, Gaussian Fisher matrix results in comparison to 
our A^-body results for the fiducial field. Finally, we conclude 
in§|8l 

2. NUMERICAL FISHER MATRIX ANALYSIS 

To study the cosmological information from the log- 
transformed field numerically, we use a large set of ray- 
tracing simulations. The ray-tracing simulations are con- 
structed from 2 X 200 realizations of A^-body simulations with 
box sizes of 240 and 480/i^^Mpc on a side, respectively. The 
number of particles in each simulation is 256^. For the fidu- 
cial cosmology, we adopt the concordance ACDM model: 
matter fraction il„i = 0.238, baryon fraction fif, — 0.042, 
dark energy fraction fix = 0.762 (therefore a flat universe), 
the equation of state parameters wq — —1 and Wa = 0, spec- 
tral index — 0.958, normalization As = 2.35 x 10~^, and 
Hubble parameter h = 0.732. Note that the fiducial cosmol- 
ogy gives (78 = 0.76 (the variance of the present-day density 
fluctuation in a sphere of radius 8 h^^Mpc). We assume three 
delta-function like source redshifts at z = 0.6, 1.0, and 1.5 
for a tomographic study. From the 400 A^-body simulations, 
we generate 1000 realizations of 5° x 5° lensing convergence 
fields (i.e., a total of 25000 square degrees) with 2048^ pix- 
els (0.15 arcmin per pixel) at each source redshift using ray- 
tracing (a total of 3000 convergence field s). Details of the ray- 
tracin g can be found in lSato et al.l (l2009h (see, also lSato et al.l 
l20Tlh . 

We resample the convergence fields in 128^ pixels (2.4 ar- 
cmin per pixel) by averaging over nearby 16 by 16 pixels as 
our fiducial case. As will be discussed in § |6] we find that 
using a larger pixel leads to a more efficient logarithmic map- 
ping. We compute the power spectra of the 3 x 1000 con- 
vergence fields assuming periodic boundary conditions be- 
fore and after logarithmic mapping. The covariance matrix 
is derived by calculating covariance between band powers at 
different wavenumber bins and at different source redshifts. 
The resulting covariance matrix represents dispersions in the 
lensing power spectra for an area of 5° x 5°. We assume 
a future, wide-field weak lensing survey of 50 00 square de- 
grees , just as in the Dark Energy Survey (DES; lAhbott et all 
120051) . by rescaling each elements of the covariance matrix by 
1/(5000/25). 



For the power spectrum and the covariance calculations, we 
need to determine the multipole bin width. An important re- 
quirement is that the dimension of the covariance matrix is 
smaller than the total number of realizations that are used for 
generating it. The dimension of the covariance matrix with 
three source redshift bins is three times bigger than that with 
a single source redshift bin. We choose Al = 200 so that go- 
ing up to /niax = 2000 requires a 30 by 30 covariance matrix. 
The dimension of the covariance matrix is then smaller than 
the 1000 realizations. Another reason for using Al = 200 
rather than a smaller bin width is to reduce the sample vari- 
ance effect in the derivatives calculation and therefore to make 
the derivatives smoother 

2.2. Log-mapping 

We use the logarithmic mapping that was introduced in Pa- 
per I. A log field is defined as: 



Kin{0) = Ko In 
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where kq is a constant with a value slightly larger than the ab- 
solute value of the minimum value of k for a given cosmology, 
source redshift, and survey pixel - this keeps the argument of 
the logarithm positive. In detail, to generate the covariance 
matrix for the fiducial cosmology, we use kq = |K,„in| +0.001 
where K,nin is the minimum pixel k of the 1000 convergence 
fields at each source redshift: Umin = -0.01348, -0.02940, 
—0.04874 at source redshifts Zg = 0.6, 1, and 1.5, respec- 
tively. While Kin reduces to the standard convergence in the 
limit of small k, the log alters k in very high or low conver- 
gence regimes. The parameter kq tunes the degree of the al- 
teration such that, the smaller kq, the more we alter the field. 
Note that various properties of the Kin field that will be dis- 
cussed in this paper show a slow, asymptotic behavior as a 
function of kq near the minimum value of k. Therefore fine- 
tuning is not necessary for kq. 

2.3. Derivatives 

For the Fisher matrix calculation, we need to first calcu- 
late numerical derivatives of the lensing power spectra with 
respect to cosmological parameters. To do this, we ran ray- 
tracing simulations of the convergence fields for cosmolog- 
ical models perturbed around the fiducial cosmology. We 
varied each of the following cosmological parameters: As, 
Us, the cold dark matter density ilji^ (the fiducial values for 
rtch^ = 0.1054), Ox, and wq by ±10%, respectively, and Wa 
by ±0.5, respectively. The Hubble parameter h, 0„i, and fl^ 
are then dependent parameters. For each of the 12 different 
cosmological models, we built 40 realizations of the conver- 
gence fields for each of the three redshift bins. Therefore we 
use a total of 3(40 x 12) = 1440 simulated convergence fields 
for calculating the derivatives. We then computed power spec- 
tra with and without the logarithmic transformation from each 
realization, averaged them, and derived the derivatives by dif- 
ferencing them. The value K^in (Eq. [TJ is derived for each 
cosmology and therefore different for different cosmological 
models. Meanwhile, calculating derivatives becomes tricky 
in the presence of shape noise for the logarithmic transforma- 
tion, as will be discussed in the next section. 



2.1. / bin width 
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2.4. Shape noise 

The intrinsic ellipticities of source galaxy shapes cause 
white noise c ontamination to the lensing power spectrum 
( lKaiseij[l998h . For the simulated convergence map, we can 
include the noise contamination by adding, to each pixel, ran- 
dom Gaussian noise with variance 



- Q 7 



(2) 



where ae is the rms of intrinsic shear per component, Ug is the 
mean number density of galaxies and ftpix is the pixel area. 
We set (Tj = 0.22 and fig ~ 30 arcmin"^ as for our fidu- 
cial values at each redshift, therefore a total of 90 galaxies 
per arcmin^ , which is much deeper th an current weak len sing 
surveys such as Dark Energy Survey ( lAbbott et alj|2005l) . In 
the presence of shape noise, the fiducial, k, field is closer to a 
Gaussian field and |Kmin| becomes larger than without shape 
noise; the log-transform becomes less efficient. We measure 
Kmin in the presence of the shape noise and use it when cal- 
culating the covariance matrix of the log-transformed field. 

For the fiducial field before log mapping, we calculate the 
derivatives from the power spectra without shape noise, be- 
cause the shape noise does not depend on cosmological pa- 
rameters. That is, the shape noise is included only in the co- 
variance matrix. This is consistent with what we would do 
in the likelihood analysis with real data. We would derive a 
signal po wer spectrum after s ubtracting the shape noise con- 
tribution (iHikage et al.ll20lTl) and compare it with a theory 
power spectrum as a function of cosmology. 

Including the shape noise effect in the derivatives calcula- 
tion is non-trivial for the log-transformed field. The shape 
noise not only increases noise in the measured convergence 
field, but also makes the log-transform less efficient by in- 
creasing That is, when a convergence field is op- 
timally log-transformed by using K,„i„, the power spectrum 
is different for a different level of shape noise even after the 
shape noise is, ideally, subtracted. In addition, while we con- 
duct a nominal subtraction of the shape noise from the power 
spectrum by subtracting a constant number that is derived 
from the difference in the variance of the field before and af- 
ter the log-transform in the presence of shape noise, in truth, 
the effect of the shape noise is neither scale-independent nor 
cosmology-independent in the log-transformed field due to 
the nonlinear transformation. Therefore we need to derive the 
derivatives differently for different cosmologies for the log- 
transformed field. For each of the 12 cosmological model, 
we first estimated the K,„in value for the log-transformed field 
with the shape noise. In order to minimize the impact of the 
random noise on the derivative calculation, for each cosmol- 
ogy, we cloned each of the 40 realizations 25 times using a 
different random seed for the shape noise. This way we gen- 
erate 1000 realizations from the underlying 40 realizations 
for each cosmology; the averaged power spectrum will have 
much smaller sample variance on the effect of shape noise. 
We then log-transform the fields and differentiate the result- 
ing power spectra between different cosmologies and calcu- 
lated the derivatives with respect to cosmological parameters. 

2.5. Fisher matrix 

We want to propagate the errors on the convergence power 
spectrum, our observable, into the projections of cosmologi- 
cal parameters using a Fisher information matrix formalism. 
Using the numerical nonlinear derivatives and the measured 



covariance matrices, we can compute the Fisher informa- 
tion matrix for cosmological param eters of our interest, pi, 
jTegmarkll 1 99T[ iTegmark et al.ll 1 9971) : 
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dpi 
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where Cz, (l) is the measured convergence power spectrum at 
multipole bin / and Zg', Cov is the measured covariance matrix 
between wavenumbers and source redshifts. The set of cos- 
mological parameters is: the amplitude of the power spectrum 
(in As), the slope of the primordial spectrum (ji,,), the physi- 
cal matter density in units of the critical density {Qjjih"), the 
dark energy density fix, and two parameters for the dark en- 
ergy equation of state: w = wq + z/{l + z)'Wa- Note that, by 
incorporating the information from the three source redshifts 
into the Fisher matrix, we are including the tom ographic lens- 
ing information ( lHull999tlTakada & Jainl2()04h . The error on 
the i-th parameter including marginalization over uncertain- 
ties in other parameters is estimated as (T{pi) ~ \/ {F~^)ii, 
where is the inverse of the Fisher matrix. The un- 

marginalized error is given as (T{pi) = l/\/Ei^. We have 
also tried including independent, free shape noise as param- 
eters for the three redshift bins, but find little effect on the 
constraints we derive. 

2.6. Planck prior 

Lensing information alone cannot determine all the cos- 
mological parameters simultaneously due to severe parameter 
degeneracies. We therefore include in these projections the 
CMB information expected from the Planck experiment. The 
CMB is sensitive to two additional parameters, the baryon 
density ftbh^ and curvature density Hk- We combi ne the con- 
vergence field information with the Planck prior (lHull2002h 
by adding the 8x8 Planck Fisher matrix to the 6x6 dimen- 
sional convergence Fisher matrix: 



F, 



■WL 



;icmb 
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where the lensing Fisher matrix F^^ has non-zero entries in 

the 6x6 block associated with In As, Ug, flch'^ (or ilmh'^), 
^Ix, Wq, and Wa, and the Planck Fisher matrix F™* has non- 
zero entries in all 8 x 8 elements including additional two pa- 
rameters: flijh^ and Hk- The Planck Fisher matrix we use in- 
cludes marginalization over uncertainties in the optical depth 
parameter r on which the CMB power spectra depend. The 
zero Hk and ftbh^ elements in the convergence Fisher ma- 
trix reflect the fact that these two parameters are constrained 
solely by the CMB information. 

3. GENERAL PROPERTIES OF THE LOG 
TRANSFORMATION 

Before presenting the Fisher matrix analysis, we first revisit 
the general properties of the log-transform that have been dis- 
cussed in Paper I, such as the power spectra, 1 -point prob- 
ability distribution (PDF), the information content, and the 
structure of the covariance matrix for the three source redshift 
bins. 

The upper left panel of Figure[T] shows the measured power 
spectra of the convergence field k (open symbols and/or 
dashed lines) and the log-transformed field Kin (solid sym- 
bols and/or soUd lines) at Zg — 0.6 (triangles and/or black). 
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Fig. 1. — Upper left: power spectra at z = 0.6 (black), 1 (blue), and 1.5 (red) before (dashed lines/open points) and after the log-transform (soUd lines/filled 
points). Upper light: PDF distiibution before and after the log-transform. Lower left: the information content of the individual redshift bins. Lower Right: the 
information content from the combination of the three Zs bins. Since the three bins share information, the combination of the three bins (solid and dashed lines 
with points) show less improvement after the logaiithmic transformation. Dotted lines show the Gaussian limit. 



1 (squares and/or blue), and 1.5 (circles and/or red). As ex- 
pected, the log-transform reduces the small-scale nonlinear 
clustering. The reduction is more prominent for the lower 
source redshift, the reason for which is apparent in the PDF 
distribution in the upper right panel: the PDF of k deviates 
more from a Gaussian PDF at Zg = 0.6 and therefore is more 
improved by the log-mapping. The lower left panel shows 
the improvement in the information content, the cumulative 
signal-to-noise ratio (S/N)'^ integrated up to a given maxi- 
mum multipole ^max, for each of the three Zg bins, defined 
as: 



(^max) 



J2 CiCov-\lJ')Ci, 



i.i'<i„ 



(5) 



where C; is the power spectrum of multipole I before and af- 
ter the log-transform. Gov is the covariance matrix describing 
correlations between the power spectra of multipoles I and /' 
(/, r < Zmax) at each Zg, and the summation runs over all the 



multipoles I and V su bject to 1,1' < l^ax (ISato et al.ll2009t 
iTakahashi et al.ll2"009l) . As pointed out in Paper I, this infor- 
mation content can be understood as the inverse of the frac- 
tional error on the amplitude of the observed, nonlinear power 
spectrum before and after the logarithmic transform. The 
Fisher matrix analysis will show to what extent this improved 
fractional error on the observed amplitude remains when a 
full set of cosmological parameters is used. In the lower left 
panel, the improvement in the information content due to the 
log-transform is largest for Zs = 0.6 and smallest for — 1.5 
due to the level of nonlinearities, as expected from the up- 
per two panels. The dotted gray line is the (S/N)'^ expected 
for a Gaussian case. Due to the nonlinear structure growth 
that causes significant off-diagonal covariances, the measured 
(S/N)'^ values of the k field are much smaller than the Gaus- 
sian limit. We find that the K\a field returns the (S/N)'^ closer 
to that of the Gaussian case, which confirms the results in Pa- 
per I but using a different A/. The improvement is a factor 
of ~ 6.9, 3.9, and 2.2 at Zg = 0.6, 1, and 1.5, respectively, 
at ^max ^ 1000 and a factor of 12, 7.9, 4.6, respectively at 
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Fig. 2. — Slices of the covariance matrix to show the covariance matrix property before (black dashed) and after (red solid) the log transform. Each panel 
shows a different block of the covariance matrix constructed by pairing two of the three Zs bins. The elements are normalized relative to the diagonal elements, 
i.e., Cov{l, I')/ (Cov{l, l)Cov{l' , r))^/^. In the notation Covij, i and j indicates a source redshift bin. For example, Covn is the normahzed auto covariance 
matrix of Zs = 0.6. For each block of the covariance matrix, we find that the Kin fisld is much more diagonal than ft. That is, the covariance between different I 
bins is significantly reduced after the log-transform, while the covariance between different Zs given I remains similar before and after the log-transform. 



^max 2000; note that the improvement is largest for the 
lowest source redshift. The improvement is slightly better 
than was reported in Paper I, which is mainly due to the loga- 
rithmic transformation seemingly being more efficient for the 
larger pixels used here. The lower right panel shows the im- 
provement in the information content when the information 
from the three Zs bins is combined. This is done by includ- 
ing the measured covariance between different Zg bins in Gov 
in Eq. ^ and summing the signal-to-noise ratios up to Zmax 
and over all the three Zg bins. Since the three Zg bins share 
some of the lensing structures, there are non-vanishing co- 
variances between the power spectra of different Zg bins and 
therefore the improvement due to the log-transform is some- 
what smaller when all Zg are combined: 1.9 at /max = 1000 
and 2.7 at /max = 2000. The Gaussian case in the lower right 
panel is derived also taking into account the expected covari- 
ance between different redshift bins. 



In addition to the improved information content, one of 
the potentially advantageous features of the log-mapping is 
the improvement in the covariance matrix property: it re- 
duces the size of the off-diagonal terms (Paper I). Figure |2] 
shows two rows of the covariance matrix for / = 306 and 
/ = 1100. Each panel shows a different block of the nor- 
malized covariance matrix constructed by pairing two of the 
three Zg bins. For each block of the covariance matrix, we 
find that the Kin field is much more diagonal than k. That 
is, the covariance between different / bins is significantly 
reduced by the log-transform. Both k and Kin fields show 
a slightly higher level of off-diagonal covariance compared 
to what we found in Paper I, partly due to the larger Al 
bin used here (i.e., Al = 200 compared to Al = 100 in 
Paper I) and also probably due to sample variance. Note 
that the different bin width alters only the Gaussian covari- 
ance contribution dMeiksin & Whitelll999l : [Scoccimarro et alj 
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Fig . 3 . — Derivatives of convergence power spectra of k and fields with respect to various cosmological parameters in absolute values. Points are derivatives 
estimated from A'-body simulations; dotted lines use Halofit (for the fiducial map) as a comparison. Left: the fiducial map. Right: the log-mapping. The fractional 
derivatives (bottom panels) remain similar even after the log-transform for some parameters, especially at high source redshifts. The fractional derivatives are 
rescaled by the factor denoted in the legends for clarity. 



[Tggi lCoorav&Hull200Tt iTakada & Jainll2009l) : the larger 
bin width reduces the Gaussian covariance, the diagonal com- 
ponents of the covariance matrix, and thus increases the rela- 
tive off-diagonal components. While the covariance between 
different wavenumbers is decreased, the covariance between 
different Zg given / remains similar before and after the log- 
transform by looking at the location of the peaks in Covy for 

In summary, we observe general properties of the log- 
transform that are consistent with the results in Paper I: the 
1 -point PDF is more Gaussian, the covariance matrix is closer 
to a diagonal matrix, and the information content is greatly 
improved after the log-transform. We next propagate this im- 
provement to the errors on cosmological parameters using the 
Fisher matrix formalism. 

4. FISHER MATRIX ANALYSIS 

We study how the improvement on the precision of the am- 
plitude, or the information content, propagates into the pre- 
cision of cosmological parameters. The improvement in the 
information content for the K\a field is due to the improved 
properties of its covariance matrix. The Fisher matrix formal- 
ism then combines this with an extra piece: the sensitivity of 
the K and Kin power spectra to cosmological parameters. 

Before presenting the Fisher matrix results, we take a look 
at the numerical derivatives calculated from the convergence 
fields from a large set of A^-body simulations. Figure[3]shows 
the derivatives of the power spectrum of k (left panels) and 



Kin (right panels) fields with respect to various cosmologi- 
cal parameters (square points). As a com parison, the dotted 
lines show the prediction based on Halofit (ISmith et al.ll2003h 
for the fiducial mapping. We see an obvious decrease in the 
relative amplitude of the derivatives due to the log-transform 
(bottom panels) at Zs = 0.6. Such changes in the derivatives 
will be combined with the changes in the covariance matrix 
property in the Fisher matrix calculation. As a caveat, we find 
the relative amplitude of the derivativ e with respect to n ., de- 
creases after the log-transform, while iNeyrinckl (120 1 1 d) finds 
an increase on small scales; it might be due to a difference 
in details of the log-transform between the weak lensing field 
and the density field. 

4.1. Without Planck prior 

We first show the Fisher matrix results of the convergence 
field without shape noise contamination and CMB informa- 
tion in Figure |4] We use all the information up to ^ = 2000. 
For reference, if the amplitude of the power spectrum were 
the only parameter (as in Paper I), the log field would lead to 
an error A\n{As) = 1.9 x 10~^ for these survey parameters 
while the standard k field would have A In(yls) = 2.9 x 10^'^. 
We call this a factor of 1.5 improvement in the ID Figure of 
Merit. When we generalize to 6 parameters, the correspond- 
ing Figure of Merit (FoM, hereafter) is the square root of the 
determinant of the 6x6 Fisher matrix. In this case, we find a 
factor of 23 improvement. This is better than the naive expec- 
tation of 1.5^ — 11, so our first conclusion is that the advan- 
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TABLE 1 

Marginalized ERRORS on each parameters. 





In As 


Tls 






wo 








K alone 


0.0642 


0.897E-02 


0.493E-02 


0.472E-02 


0.458E-01 


0.119 






Kin alone 


0.0408 


0.674E-02 


0.325E-02 


0.202E-02 


0.189E-01 


0.124 






ft + Planck 


0.0101 


0.336E-02 


0.115E-02 


0.367E-02 


0.357E-01 


0.997E-01 


0.150E-03 


0.272E-02 


Kin + Planck 


0.0102 


0.304E-02 


0.107E-02 


0.149E-02 


0.161E-01 


0.632E-01 


0.148E-03 


0.253E-02 



Note. — Marginalized errors before and after the log-transform tising power spectrum information up to /max — 
2000. Both eiTors show distinct improvement after log-transform; the improvement is mainly on Q,x , wo and Wa once 
the Planck prior is included. 




Fig. 4. — Marginalized error contours before (solid black lines) and after 
log-transform (red lines) without shape noise contamination and CMB infor- 
mation. That is, we took the 2x2 sub-Covariance matrix and plotted the error 
contours. We use power spectrum information up to i = 2000. One sees that 
the marginalized errors are overall smaller for ki „ . The dotted contours show 
constraints from the Plank alone. The Planck contours for pairs between Qx . 
Wo, and Wa do not show up here due to the extreme degeneracies. 

tages of the log estimator hold up - or even increase - when 
generalizing to multiple cosmological parameters. 

Figure|4]shows some 2D slices of these constraints. The fig- 
ure shows the marginalized 1 — a error contours from weak 
lensing tomography alone for various pairs of the 6 cosmolog- 
ical parameter^ The solid black lines show the result from 
the K field before the log mapping and the red lines show the 
results of the K\n field. The constraints are tighter for the Kin 
field: the error ellipses have shrunk and the projections are 
that the K\n field often leads to narrower allowed regions. Ta- 
ble [T] lists the marginalized errors without the Planck prior: 
the log-transform shows improvement, especially on ilx and 
wq. The dotted contours in Figure |4] show constraints from 
Planck alone: the Planck contours for pairs between ilx, w^, 
and Wa do not show up here due to the extreme degeneracies. 
For all parameters other than the dark energy parameters, the 
information from the Planck mission dominates. 
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Fig. 5. — Marginalized error contours before (solid black lines) and after 
log-transform (red lines) when we combine the weak lensing data with the 
Planck mission. We use power spectrum information up to ( = 2000. One 
sees that ki^ improves constraints mainly on the three dark energy parame- 
ters, i.e., Qx, ^0, and Wa- 



of the error contours on the 8 cosmological parameters when 
we combine the weak lensing data with Planck (left panel). 
The solid black line is for the k field and the red line for the 
K\n field after the Planck Fisher matrix is combined. The dot- 
ted contour shows the error contours for the Planck alone as 
a comparison. When the Planck priors are included, the in- 
formation on parameters other than fix, wq, and Wa is dom- 
inated by the Planck information, as evident in Figure HJ /tin 
improves constraints mainly on the three dark energy param- 
eters 0. Table [T] shows that we achieve an improvement by 
a factor of 2.5 for fix, 2.2 for wq, and 1.6 for Wa by log- 
transformation. 

4.3. Improvement in the dark energy FoM 



4.2. With Planck prior 

The situation is qualitatively similar when the Planck prior 
is added in. The Figure of Merit (now 8D) is larger by a factor 
of 12 when the log estimator is used. Figure |3 shows slices 

' I.e., we take the 2x2 sub-Covariance matrix and plot the error contours. 



* The solid contours in Figure [5] do not exactly agree with the naive error 
contour combination of the solid contours in Figure [4] and the dotted Planck 
contour: we would expect this agreement if we add the Fisher matrices from 
the 2 by 2 sub-covariance matrices of the convergence field and the Planck 
mission. However we are adding the two full Fisher matrices here such that 
a given marginalized en'or contour is affected by the effect of the Planck 
information on the rest of the parameters. 
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Fig. 6. — Figure of Merit in wo-Wa (Left) and in the 8-D (middle) before and after the log-transform. The magenta and gray lines with squares in the left panel 
show the result without the Planck prior while the red and black triangles show the values with the Plank prior. We find more than a factor of 2 improvement 
in FoM„„_„^. The middle panel shows a 8-dimensional figure-of-merit, i.e., the inverse volume of the 8-dimensional parameter space. The right panel shows 
the ratio of the 8-D FoM between with and without the log-transform. We find a factor of 7-12 improvement using (max = 1000 — 2000 in the presence of the 
Planck prior. 




100 1000 10^ 100 1000 10^ 100 1000 10^ 

L,. 1™, L,. 



Fig. 7. — Left: the information content in the presence of shape noise (fig = SOarcmin"'^ for each redshift bin). The dotted fine shows the Gaussian Hmit 
in the presence of shape noise. Middle: the information content from the combination of the three source redshift bins. Right: the 8D FoM. We find very little 
improvement in the constraints on the final cosmological parameters once shape noise is added. 



We can also quantify the improvement due to the log- 
transform using the Figure-of-Merit of dark energy parame- 
ters (hereafter 'FoMwo-w.-,') that is often used in literature to 
characterize the performance of a dark energy s urvey mission. 
The D ark Energy Task Force (DETF) FoM dAlbrecht et all 
120061) is defined as 



a{wp)a{wa) ^dct(Cov[wo,u'a])' 

where Wp is the dark energy equation of state at the "pivot" 
redshift, at which the dark energy equation of state is best con- 
strained by given observables, and Cov[u'o, Wa] is the 2x2 
sub-matrix of the inverted Fisher matrix, F~^, including only 
its elements of wo and Wa- The FoM is proportional to the 
area of the marginalized error ellipse in ujq and Wa parame- 
ter sub-space. In Figure |6] the left panel shows FoMwo-w^ 
before (dashed lines) and after the log-transform (solid lines): 
we observe approximately a factor of two improvement when 
including the Planck prior (red and black lines with triangles). 

In the middle panel, we show the 8-dimensional figure-of- 
merit as a function of Imax- We find a factor of 7-12 im- 
provement using Zniax = 1000 — 2000, which can be read out 
from the right panel where we show the ratio of FoMgc be- 
tween with and without the log-transform (black line). Table|2] 
presents the FoM values as a function of /max = 1000 — 2000. 



TABLE 2 
Figure of Merit. 



^max 


FoMwo-wa 




I'oMgD 








Kin 


K 




1000 


301 


803 


3.56E-H20 


2.67E-I-21 


2000 


610 


1423 


1.90E+21 


2.21E-I-22 



Note. — Figure of Merit in wo-Wa and in the 8-D 
before and after the log-transform. 



5. FISHER MATRIX ANALYSIS WITH SHAPE NOISE 

Any weak lensing survey contains shape noise, i.e., uncer- 
tainty associated with intrinsic shapes of the galaxies. This 
noise decreases with increasing galaxy number density. In the 
presence of large shape noise, the observed field is closer to 
Gaussian (assuming Gaussian shape noise) and Kmin is larger 
due to the additional dispersion: we therefore expect that a 
log-mapping will be less efficient for a larger shape noise. 
Paper I have shown that, with a galaxy number density of 
30 arcmin"^ at = 1 and a pixel size of 2.4 arcmin, a factor 
of the improvement is 2.4 for Z,nax = 2000; we find a sim- 
ilar result of 2.3 as shown in the left panel of Figure [T] We 
assume the total mean number density of 90 arcmin^^ and 
the number densities of 30, 30, and 30 arcmin^'^ for source 
redshifts of z,, — 0.6, 1.0, and 1.5, respectively. Note that this 
is a galaxy number density that is much higher than ground- 
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Fig. 8. — In the presence of shape noise (fig = 30arcmin~ for each 
redshift bin), we find little improvement in the constraints on the final cos- 
mological parameters after log-transform. Black: before log-transform. Red; 
after log-transform. 



based, near-future weak lensing surveys. 

Figure |7] shows that, in the presence of this modest shape 
noise, the improvement in the information content due to the 
log-mapping is still large, but the 8D FoM is very similar for 
the standard k estimator and for Kin- The improvement on the 
full set of cosmological parameters is only a factor of 1.7 for 

/max = 2000. 

Figure |8] shows the marginalized error contours with the 
Planck prior in the presence of shape noise: the improvement 
due to the log-mapping is not found here despite the large 
degree of improvement in the information content. There- 
fore, we find the improvement in the information content 
by the log-mapping does not propagate efficiently to the im- 
proved cosmological information in the presence of even an 
optimistic limit of shape noise for future weak lensing sur- 
veys. The effect of s hape noise we find is consistent with 
iJoachimi et al.l (1201 ll) despite the different set of cosmo- 
logical parameters investigated, except that they find much 
smaller improvement in the information content. In calculat- 
ing the information content (Eq. |5]l, we set the signal to be 
an averaged power spectrum of the log-transformed field af- 
ter subtracting a co nstant power as an ap proximation for the 
shape noise effect; IJoachimi et al.l (1201 Ih use a power spec- 
trum without shape noise as a signal. Since the constant power 
we subtract does not include higher order contributions that 
mingles shape noise and th e clustering signal, th e signal we 
input is higher than that of IJoachimi et al.l (1201 ll) and there- 
fore our information content is larger 

It is quite possible that the problem of the shape noise lies 
in the estimator In the absence of shape noise, the log- 
transformed field is the obvious way to make the field nearly 
Gaussian. In the presence of shape noise, it is quite possible 
that one must work harder to find an estimator that recaptures 
the information lost to higher point functions. A different 



Gaussianization scheme, particularly one not so sensitive to 
Kmin, might well work better It is also possible that the in- 
formation is hopelessly lost due to the shape noise and there 
is not much to recover with log-transform. The encouraging 
results to date in the absence of shape noise suggest that this 
is an important avenue of research. 

6. THE EFFECT OF THE PIXELIZATION 

Using a larger pixel appears more efficient in terms of the 
information content than using a smaller pixel, whether with 
and without shape noise, partly due to the smaller |Kmin| that 
we can reach for a larger pixel. The top panels of Figure |9] 
shows the information content and the 8D FoM using a pixel 
of 0.6 arcmin, instead of our fiducial pixel of 2.4 arcmin, in 
the presence of no shape noise. From the left panel, one sees 
that the information content after log-mapping (left panels) 
is much lower when using a smaller pixel t han when using a 
pixel of 2.4 arcmin (in Figure [T]) (also see lNevrinckll201 m 
jjoachimi et al.| [20TTb . On the other hand, we find that the dif- 
ference in FoMgc for different pixel sizes is not as drastic 
as we have expected based on the information content result. 
The left panel of Figure[TO]shows the marginalized error con- 
tours for the pixel of 0.6 arcmin without shape noise and with- 
out the Planck prior Since this case is less efficient for the 
log-transform, unlike the case of a 2.4 arcmin pixel (in Figure 
m, the changes in the error ellipses are not a simple shrinkage 
of ellipses. The projections are that the Kin field leads to nar- 
rower allowed regions, which implies stronger degeneracies 
after the log-transform; the errors on cosmological parameters 
do not decrease after the log-transform despite the improve- 
ment in the information content. When the Planck priors are 
combined (right), the strong degeneracies in the power spec- 
trum of Kin are lifted, and now the thin error ellipses observed 
in the left panel finally translate to smaller errors on ilx, wq, 
and Wa ■ 

The bottom panels of Figure |9] shows the information con- 
tent and the 8D FoM using a pixel of 0.6 arcmin in the pres- 
ence of shape noise. In comparison to Figure|2]for pixel of 2.4 
arcmin, we find even less improvement with the log-transform 
due to a smaller pixel. The effect of pixel we observe ap- 
pears c onsistent with the effect of smoothing in lJoachimi et alj 
(l20Tll) . 

7. ANALYTIC FISHER MATRIX RESULTS 

We compare our A^-body Fisher matrix r esults w ith the an- 
alytic Fisher matrix results using a Halofit ( iSmith e t al. 2003|) 
for the fiducial mapping. As shown in the left panel of Figure 
[3] the derivatives derived from the Halofit slightly deviates 
from the A^-body results, especially on high I, such that those 
for the dark energy parameters appear similar to what was 
shown in Casarini et al. (201 1). Due to the discrepancy in the 
nonlinear convergence power spec trum between the Halofit 
and the A^-body result (lEifleiil201lT) . there seems to be a big- 
ger difference in the fractional derivatives (bottom left panel). 
We conduct a Gaussian Fisher matrix analysis using Halofit 
results and the Gaussian assumption: we call this 'Gaussian 
Fisher matrix results' . 

We find the A^-body Fisher result predicts better constraints 
on dark energy par ameters than the Gau ssian Fisher result, 
which is contrary to lCasarini et al.l ( 1201 ll) . In detail, the non- 
linear covariance matrix increases the error bars, relative to 
the Gaussian case, which is expected; however, the A^-body 
results appear to show less degeneracies between dark energy 
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Fig . 9 . — Effects of pixel size. We use a pixel of 0.6 arcmin for this figure. Top: without shape noise. Bottom: with shape noise {rig = SOarcmin ^ for each 
pixel). Left: the information content. Right: improvements due to log-mapping in terms of various quantities. 



parameters (fix jWo,Wa) than predicted by Halofit. As a re- 
sult, the A^-body Fisher matrix analysis gives better constraints 
on dark energy parameters than the Halofit-based Gaussian 
case. If we hold Wa fixed or ilx fixed, the A^-body FoM be- 
comes worse than the Halofit result. In the presence of shape 
noise, the nonlinear covariance matrix is closer to the Gaus- 
sian one. Therefore the constraints from the A^-body result are 
better than the Halofit-based Gaussian case to a larger extent. 
Table|3]lists our results. 

We comment that sample variance introduces a noisy fea- 
ture in the A^-body derivatives. We have used = 200 rather 
than a smaller bin width to reduce such noisy feature. How- 
ever, any remaining noise might have affected the A^-body re- 
sults. We tried smoothing the derivatives, which did not de- 
crease the constrains from the A^-body results and therefore 
did not reverse our finding. 

8. CONCLUSION 

We have used the Fisher matrix formalism to test the impact 
of the log-transform on the cosmological parameters. We find 
that the log mapping performs much better than the fiducial 



mapping: when Planck mission information is included as a 
prior, the log-transformed field greatly improves constraints 
especially on dark energy parameters such as fix, wq, and 
Wa- In the presence of shape noise, however, the advantage of 
the log mapping quickly diminishes. We find little improve- 
ment on the cosmological parameters after log-transform even 
with fig = 30arcmin^^ at each of the three source redshift 
bins. We find the information content is not necessarily a 
good probe of the actual precision on the final cosmological 
parameters. We also find using a larger pixel allows a more 
efficient log-transform with and without shape noise. Finally, 
we find that, for the fiducial mapping, the Halofit-based Gaus- 
sian Fisher matrix calculation gives worse constraints on dark 
energy parameters than the full A^-body Fisher matrix result. 
This appears to be due to less degeneracies among dark en- 
ergy parameters in the A^-body power spectra. 

We greatly appreciate the extremely helpful comments from 
Bhuvnesh Jain. We thank Patrick McDonald and Jan M. 
Kratochvil for useful discussions. H-JS is supported by the 
Berkeley Center for Cosmological Physics. MS is supported 
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Fig. 10. — Effects of pixel size. Using a pixel of 0.6 arcmin. Marginalized error contours before (solid black lines) and after log-transform (red lines) without 
shape noise contamination. Left: without the Planck prior. Right: with the Planck prior We use power spectrum information up to i = 2000. One sees that the 
log-transform is less efficient with a smaller pixel. 



TABLE 3 

A'-BODY Fisher ANALYSIS vs analytic Gaussian Fisher 



ANALYSIS. 




ilx Wa Wa 


ft -1- Planck 
Analytic 

K + Planck 
Analytic 

K + Planck 
Analytic 


0.367E-02 0.357E-01 0.997E-01 

0.372E-02 0.451E-01 0.128 

0.335E-02 0.165E-01 

0.154E-02 0.154E-01 

0.294E-01 0.909E-O1 
0.171E-01 0.530E-01 


re -1- Planck, shape noise 
Analytic, shape noise 

re -1- Planck, shape noise 
Analytic, shape noise 

re -1- Planck, shape noise 
Analytic, shape noise 


0.626E-02 0.808E-01 0.220 

0.134E-01 0.169 0.409 

0.525E-02 0.211E-01 

0.429E-02 0.237E-01 

0.603E-01 0.185 
0.409E-01 0.131 







Note. — The analytic Gaussian Fisher results are derived us- 
ing the Halofit results and the Gaussian assumption. Shape noise 
is based on ng = 30 arcmin^^ for each source redshift bin. The 
results on parameters other than dark energy parameters are very 
similar between the two methods, mainly because these are domi- 
nated by the Planck information; we do not show these parameters 
in this table for simplicity. The rows with "-" means that the cor- 
responding parameter is held fixed. 
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